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Abstract 


The  equations  of  motion  for  a  cantilever  beam  in  bending 
vibration  are  developed  in  state  vector  form  using  a  normal 
mode  approximation.  A  linear  optimal  control  system  generates 
a  feedback  control  proportional  to  the  state  which  is  repre¬ 
sented  by  modal  amplitudes  and  velocities  determined  using 
position  information  from  sensors.  The  observer  gain  matrix 
and  the  feedback  control  gain  matrix  are  both  determined  from 
a  steady  state  optimal  regulator  which  minimizes  the  related 
quadratic  performance  index.  Control  is  applied  through 
point  force  actuators.  The  effects  on  beam  response  of  sensor 
and  actuator  location,  reduced  observation  and  control  spill¬ 
over,  restructured  control  performance  index,  and  a  reduced 
order  optimal  regulator  are  examined.  Parameter  variations 
of  the  modal  amplitudes  were  tested  to  note  effects  on  system 
stability  with  both  stable  and  unstable  controllers.  Singular 
perturbation  theory  is  used  to  reconstruct  the  control  per¬ 
formance  index,  essentially  adding  a  penalty  function  against 
any  control  vector  that  lies  in  the  suppressed  modes  subspace. 
Singular  value  decomposition  is  used  to  construct  a  transfor¬ 
mation  matrix  through  which  observation  and/or  control  spill¬ 
over  may  be  eliminated.  A  means  by  which  optimal  actuator 
locations  may  be  determined  is  also  presented. 

System  response  is  shown  to  be  very  sensitive  to  actuator 


vi 


location.  Singular  perturbation  provided  a  method  through 
which  control  spillover  could  be  minimized  for  the  actuator 
locations  chosen,  however,  it  did  not  provide  the  means  by 
which  the  actuators  could  be  positioned  such  that  the  spill¬ 
over  effect  could  be  eliminated.  Optimal  actuator  locations 
could  be  found  using  singular  value  decomposition.  Control 
and  observation  spillover  effects  could  then  be  completely 
eliminated  only  where  the  number  of  actuators  was  equal  to  or 
greater  than  the  number  of  modes  to  be  suppressed.  Robustness 
of  system  response  to  modal  amplitude  errors  was  very  good, 
and  the  unstable  controllers  tested  did  not  appear  to  seriously 
affect  this  robustness. 
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MODERN  OPTIMAL  CONTROL  METHODS 
APPLIED  IN  ACTIVE  CONTROL  OF  A 
CANTILEVER  BEAM  IN  BENDING  VIBRATION 

Introduction 

The  control  of  large  space  structures  has  received  much 
attention  since  the  advent  of  the  space  shuttle  program.  Ac¬ 
tive  control  of  the  structural  vibrations  associated  with  large 
structures  using  modern  optimal  control  methods  is  one  area 
presently  being  investigated.  The  advantages  of  using  these 
methods  lie  in  the  fact  that  the  design  process  can  be  fully 
automated  and  adapted  to  a  large  class  of  problems  rather 
than  a  specific  model  configuration.  An  optimal  number  and 
placement  of  sensors  and  actuators  could  also  be  determined. 
Despite  these  and  other  advantages,  these  methods  have  seen 
little  application  in  the  control  of  large  space  structures 
due  to  inherent  problem  areas  associated  with  using  them. 

Two  of  the  main  problem  areas  deal  with  model  reduction  and 
modelling  inaccuracies.  A  very  large  number  of  modes  may  be 
required  to  model  a  flexible  structure.  Active  control  of 
all  these  modes  would  be  impractical  due  to  computer,  sensor, 
and  actuator  requirements,  therefore  a  reduced  order  model 
of  the  system  must  be  considered,  giving  rise  to  model  reduc¬ 
tion  error.  In  a  complex  structure , modelling  can  be  accom¬ 
plished  using  finite  element  approximations,  however,  this 
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approximation  technique  produces  modal  data  whose  accuracy 
is  uncertain.  Modelling  inaccuracies  cause  the  modal  fre¬ 
quencies  and  mode  shapes  to  be  in  error. 

Balas  (Ref  1)  discusses  a  method  by  which  a  class  of 
flexible  structures  can  be  controlled  through  the  use  of 
sensors,  actuators,  a  state  estimator,  and  state  variable 
feedback;  and  he  showed  that  this  could  lead  to  closed  loop 
instabilities  in  the  reduced  model  when  modes  other  th; 
those  controlled  are  considered.  These  instabilities  are 
caused  by  the  combined  effect  of  what  Balas  termed  observa¬ 
tion  and  control  spillover.  Sesak  (Ref  2)  used  singular 
perturbation  theory  and  a  highly  mathematical  approach  toward 
eliminating  control  spillover.  Sesak' s  technique,  however, 
gave  little  insight  into  how  or  when  it  might  work.  Several 
specific  examples  were  shown,  but  no  general  results  were 
obtained.  In  an  attempt  to  explain  the  singular  perturbation 
approach,  Coradetti  presented  an  analysis  of  the  limiting 
case  of  the  singular  perturbation  method  and  showed  that  it 
was  equivalent  to  finding  a  transformation  matrix  which  re¬ 
duced  control  spillover.  This  transformation  matrix  can  be 
determined  from  singular  value  decomposition  of  the  control 
matrix  associated  with  those  modes  in  the  model  that  are  not 
actively  controlled.  Using  this  technique  actuator  locations 
could  be  predetermined  to  produce  a  transformation  matrix 
which  would  eliminate  control  spillover.  Actuator  placement 
and  the  number  of  actuators  necessary  to  achieve  this  result, 
however,  were  not  examined  in  depth;  and  applying  this 
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technique  to  eliminate  observation  spillover  was  not 
considered. 

The  purpose  of  this  thesis  is  to  apply  Balas'  method  to 
a  cantilever  beam  and  determine  the  effect  on  the  beam  re¬ 
sponse  of  sensor/actuator  location,  and  reduced  observation 
and  control  spillover  using  singular  perturbation  and  singu¬ 
lar  value  decomposition  techniques.  In  addition,  recent  work 
by  Johnson  (Ref  5)  has  shown  that  standard  optimal  control 
design  can  .lead  to  unstable  controllers.  The  effects  of  a 
stable  versus  an  unstable  controller  on  system  robustness  to 
modelling  inaccuracies,  specifically  modal  error,  is  also 
considered. 

The  approach  taken  involves  forming  a  discrete  system 
through  normal  modes  approximation.  Position  sensors  are 
used  to  determine  modal  amplitudes,  from  which  a  state  esti¬ 
mator  reconstructs  the  mode  shapes.  Control  is  accomplished 
using  point  force  actuators  through  state  variable  feedback. 
Both  the  observer  and  feedback  control  use  gains  produced 
from  steady  state  optimal  regulators.  Control  performance 
index  restructuring,  using  singular  perturbation,  penalizes 
the  control  against  acting  on  modes  other  than  those  desired. 
Singular  value  decomposition  of  the  modal  amplitudes  at  the 
sensor  or  actuator  location  is  used  to  develope  a  transforma¬ 
tion  matrix  through  which  observation  and/or  control  spill¬ 
over  may  be  eliminated. 

The  main  areas  examined  in  this  report  are:  system  sensi¬ 
tivity  to  sensor  and  actuator  placement;  reducing  observation 
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or  control  spillover  using  singular  perturbation  and  singular 
value  decomposition  techniques;  the  number  of  actuators  re¬ 
quired  and  placement  necessary  to  eliminate  control  spill¬ 
over;  system  robustness  to  modal  error  for  both  stable  and 
unstable  controllers. 
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System  Model 


Equations  of  Motion 

The  continuous  differential  equation  of  motion  for  a 
cantilever  beam  in  bending  vibration  is  defined  in  Meirovitch 
(Ref  4:208,209)  as 


El  (x) 


9 2 Y  (x, t) 


9  Xu  _L  /  -  9 

9x2  9x2 


+  f  (x, t)  =  m(x) 


9  2  Y  (x,  t) 


with  boundary  conditions 


Y  (0 ,  t)  =  0 


9Y (x , t)  . 

9x  x=0 


EI(x) 


92Y (x,t) 


9x  3x2  x=L 

where  x  is  the  variable  along  the  length  of  the  beam,  Y  is 
the  vertical  deflection,  El (x)  is  the  bending  stiffness,  f 
is  the  applied  force,  L  is  the  length  of  the  beam,  and  m  is 
the  mass/unit  length.  The  case  considered  is  that  of  a 
uniform  beam  in  free  vibration  so  that 

El (x)  =  constant 

f(x,t)  =0  C 

m(x)  =  constant 

A  discrete  form  of  Eq  (1)  can  be  obtained  through  the 


substitution 


Y(x,t) 


(4) 


6 


(9) 


cosB^L  coshB^L  =  -1. 

The  mode  shapes,  are  normalized  with  respect  to  the  con¬ 

stant  mass  density,  that  is 

L 

/  m<ji  •  (x)dx  =  6^  (10) 

0  J  J 

where  6. .  is  the  Kronecker  delta.  These  mode  shapes  and  modal 
il 

amplitudes  are  then  used  to  form  a  linear  system  model  to 
which  modern  control  theory  is  applied. 


Linear  System  Model 

The  number  of  modes,  n,  in  a  structural  model  may  still 
be  very  large.  For  practical  purposes  the  controller  should 
be  concerned  with  as  few  modes  as  possible  and  still  maintain 
a  stable  system.  A  determination  of  which  modes  need  to  be 
controlled  will  be  discussed  later,  however,  for  the  time  be¬ 
ing  it  is  accepted  that  the  formulation  of  a  linear  model 
should  allow  for  this.  If  the  model  to  be  developed  is  con¬ 
sidered  to  consist  of  three  parts,  a  controlled,  a  suppressed, 
and  an  un-modelled  part;  Eq  (4)  can  be  partioned. 

Y  (x,  t)  =  Yc(x,t)  +  Ys(x,t)  +  Ymn^t)  (11) 

Where  the  suppressed  modes  are  those  over  which  control  is 
unnecessary,  however,  their  behavior  must  be  considered.  The 
un-modelled  modes  are  generally  the  higher  frequency  modes 
where  natural  damping  is  condisered  adequate  to  preclude  the 
development  of  any  instabilities. 

The  controlled  part  of  Eq  (11)  is 


Yc(x,t) 


C 

I  (x)Ui(t) 

i=l  1 


(12) 
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the  suppressed  part 


c+s 

Y  (x,t)  =  Z  4>.(x)U.(t)  (13) 

S  j=c+l  ^  ^ 

and  the  un-modelled  part 

n 

Yum(x't}  =  Z  4>k(x)Uk(t) 
k=c+s+l 

where  c  represents  the  number  of  controlled  modes,  s  is  the 
number  of  suppressed  modes,  and  n  is  the  total  number  of  modes 
in  the  model.  Further  consideration  of  the  un-modelled  modes, 
however,  is  necessary  in  that  the  controller  that  will  be  de¬ 
veloped  in  the  linear  model  will  have  no  knowledge  of  these 
modes.  The  purpose  of  discussing  the  un-modelled  modes  is  to 
demonstrate  that  the  model  reduction  process  actually  involves 
two  truncations,  the  first  to  a  finite  number  of  modes  to  des¬ 
cribe  a  structure,  and  the  second  limiting  concern  only  to 
those  considered  critical.  Both  of  these  truncations  lead  to 
inaccuracies . 


The  amplitudes  of  the  controlled  and  suppressed  modes  to¬ 
gether  with  their  rates  of  change  form  the  states 


Vs(t) 


(t)  =  ;  uiT(t)JT 

(t)  =  [ujT(t)  :  UjT(t)J 


substituting  this  into  the  discrete  form  of  the  equation  of 
motion,  the  system  can  be  modelled  by 


V  (t)  *  AV  (t)  +  B_f  (t) 
c  — c  c  — c 


Vs(t)  =  AsVs(t)  +  BsI(t) 


where  a  dash  over  a  variable  represents  a  vector,  and  a  dash 


under  a  matrix.  f(t)  is  the  vector  of  control  inputs. 
The  system  parameter  matrices  are 

0  I 


*c  = 


“Sc  -2^i“i 


(18) 


As  = 


-ft  -2£.o). 

-S  3 


(19) 


B  = 


0 

B 


(20) 


B 

— s 


0 

B, 


(21) 


where  uk,  uk  are  the  natural  frequencies  of  the  controlled  and 
suppressed  modes  respectively,  the  ft  matrices  are  diagonal  ma¬ 
trices  whose  elements  are  the  squares  of  the  natural  frequencies, 
£  is  the  damping  factor,  £  and  I  are  the  null  and  identity  ma¬ 
trices,  Bc  and  Bg  are  the  matrices  whose  columns  are  the  mode 
shapes  evaluated  at  the  actuator  locations,  Xj...x  . 


Sc  = 


5s  = 


4>  (x  )  •  •  (x  ) 
i.i  i  . a 


(x  )•••<}>  (x  ) 
c  i  a 


*c+l(x1,,,#*c+l<xa> 


+c+s(x,,-”*c+s(xa> 


(22) 


(23) 


The  sensor  output  is  given  by 


9 


\  The  null  £  portion  of  C  and  C  matrices  represent  the  rates 

I'  c  s 

of  change  of  the  mode  shapes  at  sensor  locations  and  are  0 

j  since  displacement  only  sensors  are  used.  Although  this  model 

1 

|  has  been  developed  for  a  relatively  simple  structure  using 

I  exact  mode  shapes,  the  method  could  easily  be  used  on  complex 

I 

1  structures  using  finite  element  approximations,  where  no 

t  '  f 

S  simple  continuum  description  is  available. 

Since  we  wish  to  use  state  variable  feedback,  complete 
knowledge  of  Vc(t)  is  required,  however,  the  only  knowledge 
j  available  is  that  contained  in  the  sensor  outputs.  There- 

j  fore,  a  state  estimator  is  required  in  order  to  implement 

j  the  control  system. 
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Observer  Model  and 


Feedback  Control  System 

The  state  estimator  has  the  form 

$  =  AO  (t)  +  B_f(t)  +  K  Y(t)  -  Y(t)  (29) 

c  — C  C  .  — — 

A  __ 

where  Vc(t)  is  the  estimate  of  Vc(t)  such  that 

Vc(t)  =  Vc(t)  +  e(t)  (30) 

where  e(t)  is  the  error,  the  sensor  output  is 

Y(t)  =  C^Vjt)  +  CsVg(t)  (31) 

The  output  corresponding  to  the  estimated  state  is 

Y(t)  =  C  V  (t)  (32) 

— c  c 

Using  equations  (29),  (30),  (31),  (32)  and  (16)  it  can  be 
shown  that 

e(t)  -  (£<,  -  K  Cc)e(t)  +  K  CsVs(t)  (33) 

The  observer  gain  matrix,  K,  must  be  developed  such  that  this 

estimator  error  decays  more  rapidly  than  the  system  dynamics 

for  proper  system  response.  The  eigenvalues  of  (A  -  K  C  ) 

determine  the  decay  rate.  Since  these  eigenvalues  are  the 

T  T  T 

same  as  the  eigenvalues  of  (A^  -  K  ) ,  the  equations  of 

motion  for  the  sensing  system  can  be  modelled  in  the  form 

w(t)  =  A  Tw(t)  -  C  Tg(t)  (34) 

where  ' 

g(t)  =  KTw(t)  (35) 

The  observer  gain  matrix,  K,  is  then  calculated  using  a  steady 

state  optimal  regulator,  where 

00 

J  =  h  f  w^Qgw  +  g^R^g  dt  (36) 
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(37) 


is  minimized.  The  optimal  solution  is  given  by 

KT  =  -R  ,_1C  P 
-  -nob  -s- 

where  P  is  the  solution  to  the  algebraic  matrix  Riccati 


equation 


Z  £  +  ScE  £T«ob'‘£  E  +  as  -  0 


where  R^  and  Qg  are  weighting  matrices. 

Similarly,  the  control  feedback  gain  matrix,  G,  is 
calculated  using  the  steady  state  optimal  regulator 

00  rn 

J  =  h  f  (V/FV,  +  f  Rf)dt  C 

0  c  -  C 

where  F  and  R  are  weighting  matrices,  and 

f(t)  =  GVc(t) 

The  optimal  solution  is  given  by 

G  =  -R“ 1 B  TS  (A 

”C 

where  S  solves  the  matrix  Riccati  equation 

S  A  +  A  TS  -  S  B  R-1B  TS  +  F  =  0  « 

—  — c  — c  —  —  — c—  — c  —  —  — 

using  these  results  the  equations  of  motion  for  the  system 


are  now 


V  (t)  =  (A  +  B  G)  V  (t)  +  B  G  e  (t) 
c  — c  — c —  c  ■ — c — 


Vs(t)  =  -sVs(t)  +  -s-  vc(t)  +  e(t)  (43) 

By  defining  a  composite  system  state  vector 

Z(t)  =  [vcT(t)  :  e*1*  (t)  :  VgT(t)]  (44) 

a  closed  loop  system  model  with  state  variable  feedback  con¬ 
taining  the  effects  contributed  by  suppressed  modes  can  be 
written  as 
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r 


Z(t)  = 

The  terms  K  C  and  B  G  are  the  observation  and  control  spill- 
s  s 

over  terms,  respectively.  It  is  evident  from  this  equation 
that  the  conbination  of  the  two  could  produce  instabilities 
in  the  system. 

Linear  Control  System  Represented 
In  Block  Diagram  Form 


The  equations  now 

governing 

the  system  are 

ir  =  A  V  +  B  f 

state  equation 

(46) 

Y  =  C  V 

displacement  equation 

(47) 

?  =  G  Vc 

control  equation 

(48) 

Vc  =  ACVC  +  Bj 

+  K(Y-Y) 

estimator  equation 

(49) 

Putting  the  system  into  block  diagram  form  and  separating  the 
plant  from  the  controller  yields  Figure  1  .  Johnson,  (Ref  7) 
showed  that  through  manipulation,  this  diagram  can  take  the 
form  of  Figure  2 /  from  which  it  is  easily  seen  that  the  trans¬ 
fer  function  for  the  controller  is 

r(s)  =  K (SI  -  A  -  B  G  +  KC  )_1G  (50) 

y  C  C  C 

If  any  of  the  eigenvalues  of  (Ac  +  BcG  -  KCc)  are  positive, 
the  controller  would  be  unstable.  It  is  interesting  to  note 
that  in  a  standard  optimal  design  process  an  unstable  con¬ 
troller  could  unknowingly  be  developed,  which  when  connected  to 


A  +B  G 
~~C  —  c— 


B  G 
— s— 


B  G 
— c— 

A  -KC 


B  G 
— s— 


KC 


Z(t) 


(45) 
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the  plant  would  still  produce  a  stable  system.  In  addition 
to  the  devastating  effects  this  unstable  controller  could 
have  on  a  system  if  it  were  ever  disconnected  from  the  plant, 
possibly  through  an  electrical  interuption,  a  sensitivity 
study  of  system  response  to  parameter  variation  was  conducted 
for  both  stable  and  unstable  controllers. 
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Plant 


Fig  2.  System  Represented  in  Modified  Block  Diagram  Form 
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Singular  Perturbation  Optimal  Control 

This  section  outlines  the  procedure  developed,  by  Sesak, 
of  applying  singular  perturbation  theory  to  the  control  of  a 
flexible  structure.  The  method  is  designed  to  decouple  the 
suppressed  and  the  controlled  systems  by  reducing  the  signal 
level  in  the  suppressed  portion  of  the  plant.  Under  these 
restrictions  the  controller  does  not  excite  the  suppressed 
modes,  therefore,  the  suppressed  forced  response  (control 
spillover)  goes  to  zero.  It  can  be  seen  in  Eq  (45)  that  if 
the  control  spillover  is  zero,  the  eigenvalues  of  the  system 
(  would  be  the  eigenvalues  of  the  diagonal  elements.  Since 

these  are  designed  to  be  negative  the  system  would  then  be 
stable. 

Singular  perturbation  is  a  variation  that  changes  the 
order  of  the  differential  equation  characterizing  a  particular 
dynamic  system.  Consider  the  performance  index  of  the  system 
under  consideration  given  as 

J  =  h  f  (VcTFVc  +  VST  Qs  Vs  +  fTRQf)dT  (51) 

subject  to 


It  is  desired  that  the  dynamics  associated  with  the  suppressed 
modes  be  stable,  therefore,  the  singular  perturbation  constraint 
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V  =  0  is  used.  This  yields 
s 


or 


0  =  A  V  +  B  f 
s  s  s 


V  —  -A  _1B  f 
s  as 


substituting  this  into  Eq  '(51)  yields 

VeTQeVe  =  fTBcT(Ae-1)TQs(As"1)BRf 


s  s  s 


(53) 


(54) 


(55) 


where  Qg  is  a  positive  semi-definite  weighting  matrix.  The 
new  performance  index  is 
00 

J  =  h  f  V0TFVr  +  f T ( Rn  +  B„T (Ac" 1 ) TQqAq- 1 Bs) f dt  (56) 


subject  to 


VC  =  ACVC  +  Bcf 


(57) 


The  weighting  matrix,  Qs,  is  designed  to  change  the  weighting 
of  the  penalty  on  the  suppressed  modes  until  control  spill¬ 
over  is  eliminated.  Whether  or  not  this  is  always  possible, 
or  under  what  conditions  it  is,  was  not  discussed.  Several 
specific  examples  were  shown,  but  no  general  results  were 
obtained.  It  was  found  that  the  method  would  reduce  control 
spillover,  using  very  large  gains,  however,  the  amount  of  re¬ 
duction  was  very  dependent  on  actuator  locations. 
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Transformation  Method 


A  method,  presented  by  Coradetti,  for  eliminating  con¬ 
trol  spillover,  using  a  transformation  matrix,  is  developed 
in  this  section.  The  transformation  matrix,  which  is  deter¬ 
mined  from  singular  value  decomposition  of  the  suppressed 
modal  amplitudes,  is  used  to  constrain  such  that 


and 


BSKC  =  0 
?  0 


(58) 

(59) 


thereby  eliminating  control  spillover,  and  allowing  complete 
control  over  the  modes  chosen.  The  matrix  Bs  may  be  express¬ 
ed  through  singular  value  decomposition  in  the  form 

Bg  =  U  W  XT  (60) 

where 

(m  x  n)  is  a  rectangular  matrix  of  modal  amplitudes 
evaluated  at  the  actuator  locations 
U  (m  x  m)  is  the  orthogonal  matrix  of  left  singular 
values 

X  (n  x  n)  is  the  orthogonal  matrix  of  right  singular 
values 

,  where  all  elements  are  zero  except 

°ii'  *  =  These  are  the 

non-zero  singular  values  of  Bg. 

The  number  q  is  equal  to  the  rank 

of  B  . 

— s 


W  (m  x  n)  = 


a 

0 
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(61) 

(62) 

(63) 

(64) 

(65) 


BgT  =  UqoXqTX  =  U^a  0=0  (66) 

In  order  to  acheive  Eq  (66)  certain  conditions  have  to  be  met. 

If  the  rank  of  Bg (m  x  n)  is  r,  then  the  nullspace  formed  of 
solutions  to  BgT  =  0,  has  n  -  r  free  variables  as  independent 
parameters.  If  r  =  n,  there  are  no  free  variables  and  the  null- 
space  only  contains  T  =  ().  In  the  case  of  r  <  m  there  are  n  - 
r  constraints  on  T  in  order  for  BST  =  0  to  be  solvable.  If 
one  solution  exists,  then  every  other  solution  differs  from  it 
by  a  vector  in  the  nullspace  of  B  .  Therefore,  in  the  case  of 
r  =  n  <  m  the  only  T  matrix  that  is  possible  for  BgT  =  0  is 
T  =  0.  However,  when  r  <  m,  n,  then  there  will  always  be  a  solu¬ 
tion  to  B  T  =  0,  for  which  T  ^  0.  If  B  is  of  full  rank,  r  = 
s  —  — s 

min  (m,n) ,  and  must  therefore  be  reduced  to  realize  Eq  (66) . 

Since  m  is  determined  by  the  number  of  modes  suppressed  and 


n  by  the  number  of  actuators  in  the  system,  a  T  matrix  can  be 
found  which  satisfies  Eq  (66)  only  when  the  number  of  actuators 
are  equal  to  or  greater  than  the  number  of  modes  to  be  suppressed. 
The  minimum  number  of  actuators,  that  can  be  used  successfully 
with  this  method,  is  two;  since  a  matrix  of  rank  1  cannot  be 
reduced. 

The  reduced  order  optimal  regulator  problem  is  defined  by 


J  =  %  /  votfvc 


subject  to 
with  the  additional  constraint 


vc  =  £cVc 


+  fTRf  dT 
+  Bcf 


Bsf  =  0 

Defining  a  new  control,  U,  where 

I  =  TU 

The  new  optimal  regulator  problem  is 


subject  to 
where 


J  =  h  f  Vf,TFVc  +  UTRJLJ  dT 
0  "T 

V  =  A  V  +  B  U 
c  — c  c  t 


and 


The  solution  is  given  by 


B  =  B  T 
— T  ~c~ 


Rj,  =  TTRT 


u  - 

=  Ws 

where  S  is  solved  from  the  matrix  Ricatti  equation 

0  =  F  +  S  A  +  A  TS  +  S  B  R  ']B Js 
—  —  —  — c  — c  —  —  — T—T  ~ T  — 


(67) 

(68) 

(69) 

(70) 


(71) 

(72) 

(73) 

(74) 

(75) 

(76) 

(77) 
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The  gain  matrix  is  then 


Kc  =  TKt  (78) 

Actuator  locations  were  determined  using  the  singularity 
computer  program  (Appendix  B) .  The  program  was  used  to  vary 
actuator  locations  and  determine  the  nonzero  singularities 
of  Bs.  Actuator  locations  were  chosen  such  that  one  of  these 
singularities  approached  zero,  indicating  a  reduction  of  the 
rank  in  Bs.  These  actuator  locations  were  then  used  in  the 
main  program  and  the  T  matrix  constructed  through  singular 
value  decomposition  of  Bg  at  these  locations. 


Limit  Method 


This  section  contains  an  extension  of  Coradetti's 
Transformation  Method  which  further  explains  Singular  Per¬ 
turbation  Optimal  Control.  This  information  was  not  avail¬ 
able  durring  the  investigation  conducted  on  Singular  Pertur¬ 
bation  Optimal  Control  and  it  serves  to  substantiate  the 
results  obtained. 


Considering  Eqs  (75)  and  (76)  ,  the  matrix  Ricatti  Eq 
(77)  can  be  rewritten  as 


0  =  F  +  S  A  + 
—  — - c 

A^S  +  s  bct(ttr  t)-1ttbcts 

(79) 

Defining  a  new  matrix. 

H,  where 

H  =  T(TTR  T)-1TT 

(80) 

Eq  (79)  becomes 

0  =  F  + 

S  A  +  A^TS  +  S  B  H  B 
—  — c  — c  —  —  — c—  — c  — 

(81) 

The  associated  performance  index  is 


J 


H  f 


—  T  — 
V  FV 
c  —  c 


+  fTH~1f 


dT 


(82) 


subject  to 


V  =  A  V  +  B  f  (83) 

c  — c  c  — c 

Comparing  Eqs  (81)  and  (82)  to  Eqs  (41)  and  (67),  it  can  be  seen 

that  H~ 1  =  R  and  R~ 1  =  H ;  however,  by  design,  T(TTR  T)TT  is  not 

full  rank.  Therefore,  H"1  does  not  exist;  consequently,  there 

is  no  finite  matrix  R  =  H_1  such  that  R"1  =  H,  except  possibly 

in  a  limiting  case.  Defining  a  new  control  weighting  matrix,  Rg 

R  =  R  +  aB  TQ  B  =  H-1  (84) 

— s  —  — s  —  — s  — 
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where  Q  is  any  positive  definite  matrix  and  a  is  any  positive 

real  scalar.  Coradetti  shows,  in  an  involved  mathematical 

proof,  that,  in  the  limiting  case  of  a  -*■  »,  R  -1  =  H.  Adding 

s 

this  to  the  performance  index  yields 

oo 

J  =  h  f  v  ?FV„  +  T(R  +  aBeTQ  B  )F  dT  (85) 

Q  c  -  c  -  -s - s 

It  is  now  quite  clear  that  it  is  penalizing  the  performance 
index  against  any  control  vector  that  lies  in  the  subspace 
occupied  by  B  .  Multiplying  B  TQ  B  by  a  large  enough  scalar 
would  effectively  eliminate  any  control  acting  on  the  suppress¬ 
ed  modes,  however,  in  theory,  it  would  require  an  infinite 
value  to  accomplish  this.  Coradetti  found  that  this  method 
will  decrease  control  spillover  until  the  numbers  involved 
become  so  large  that  computer  roundoff  errors  cause  problems 
in  matrix  inversion,  necessary  in  solving  the  matrix  Ricatti 
equation.  Singular  Perturbation  Optimal  Control  essentially 
applies  the  Limit  Method  with  the  substitution 

aQ  =  1 ) TQSAS" 1  (86) 

made  in  Eq  (84)  .  The  Limit  Method,  as  presented,  shows  that 
while  Singular  Perturbation  Optimal  control  can  reduce  control 
spillover,  it  can  never  actually  eliminate  it. 
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Computer  Model 

The  function  of  the  main  computer  program  was  to  assemble 
the  system  matrix  and  determine  the  associated  eigenvalues. 

The  program  was  modified  as  needed  to  incorporate  singular 
perturbation  optimal  control,  the  transformation  method,  and 
parameter  variations.  A  listing  of  the  program  using  the 
transformation  Method,  and  permitting  parameter  variation  is 
listed  in  Appendix  A. 

The  first  function  of  the  program  was  to  assemble  the 
required  matrices.  Input  data  necessary  consisted  of:  beam 
length,  mass,  and  bending  stiffness;  sensor  and  actuator  lo¬ 
cations;  and  the  roots  to  the  characteristic  equation.  The 
root  values  input  were:  mode  1  -  1.875,  mode  2  -  4.694,  mode 
3  -  7.855,  mode  4  -  10.996,  mode  5  -  14.137,  mode  6  -  17.278. 
The  beam  length,  mass,  and  bending  stiffness  were  1,  and  F, 

,  5  and  Rq].,  weighting  matrices  were  each  the  identity  ma¬ 
trix.  The  natural  frequencies,  determined  by  squaring  the 
roots  to  the  characteristic  equation,  are  then  used  to  con¬ 
struct  the  A-  and  A„  matrices.  Equation  (8) ,  for  the  modal 
amplitudes,  is  solved  using  subroutine  Modes  (Appendix  A) , 
and  used  to  fill  matrices  Bc,  Bg,  and  Cg.  Using  F,  B^, 

A c  and  R,  the  steady  state  algebraic  Ricatti  equation  is 
solved  using  subroutine  MRIC  from  the  Aerospace  Medical  Re¬ 
search  Laboratory  Library.  The  optimal  control  gain  matrix. 
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G,  is  then  solved  for  using  Eq  (40)  .  The  optimal  observer 
gain  matrix,  K,  is  solved  for  in  an  identical  manner.  The 
system  matrix,  Eq  (45) ,  was  then  constructed,  and  subroutine 
EIGRF  from  the  International  Mathematical  and  Statistical 
Library  (IMSL)  produced  the  complex  eigenvalues  associated 
with  the  system.  The  weighting  Matrix,  Qg,  was  varied  from 
1  to  1  x  10 14  in  examination  of  the  singular  perturbation 
optimal  control  method.  Subroutine  LSVDF ,  from  the  IMSL  li¬ 
brary  uses  singular  value  decomposition  to  construct  the 
transformation  matrix,  T,  from  which  a  reduced  order  gain 
was  developed. 

A  second  computer  program,  listed  in  Appendix  B,  was  used 
to  generate  data  for  determining  sensor  and  actuator  placements 
In  this  program  sensor  or  actuator  location  was  varied  and  the 
non-zero  singularities  associated  with  C  or  B  were  determined 
through  singular  value  decomposition.  The  desired  sensor  or 
actuator  locations  were  those  which  caused  one  of  the  non-zero 
singularities  to  most  closely  approximate  zero. 

The  sensitivity  study  was  made  using  up  to  +20%  variation 
of  the  elements  of  the  B  .  B  .  C  ,  and  C„  matrices.  This  was 
accomplished  by  multiplying  the  elements  of  these  matrices  by 
one  factor,  if  the  sum  of  its  subscripts  were  even,  and  by 
its  negative  counterpart  if  the  sum  was  odd.  This  effectively 
alters  the  mode  shapes,  simulating  modal  error.  Although  this 
method  does  not  account  for  any  frequency  errors  incurred,  the 
results  should  essentially  reflect  the  effects  produced  by 
modal  amplitude  error. 
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Outline  of  Investigation 

An  investigation  was  conducted  examining  Singular  Per¬ 
turbation  Optimal  Control  and  the  Transformation  Method 
applied  to  the-  model  developed  of  a  cantilever  beam.  In 
addition,  an  investigation  to  determine  system  robustness 
to  parameter  variations  was  made. 

The  first  area  investigated  was  the  use  of  Singular 
Perturbation  Optimal  Control  to  eliminate  control  spillover. 

A  single  sensor  and  actuator  were  initially  used  to  control 
three  modes  and  suppress  one.  When  it  was  found,  that  in  the 
process  of  reducing  the  control  spillover,  control  over  the 
first  three  modes  was  gradually  being  lost,  another  sensor 
and  actuator  were  added  to  the  system.  It  was  then  possible 
to  decrease  control  spillover  and  still  maintain  control.  The 
conditions  which  determine  the  reduction  realized  using  this 
method  are  discussed. 

The  second  area  investigated  was  the  Transformation 
Method  for  eliminating  control  spillover.  Two  sensors  and 
two  actuators  were  used  to  control  two  modes  and  suppress  two. 
Demand  on  the  system  was  then  increased  to  control  three  modes 
and  suppress  three.  The  effectiveness  of  the  method,  in  elimi¬ 
nating  spillover  effect  and  uncoupling  system  eigenvalues,  under 
these  conditions,  was  examined  and  some  general  conclusions 
were  drawn. 
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The  elimination  of  observation  spillover  using  the 
Transformation  Method  was  then  investigated. 

Finally,  an  investigation  of  system  robustness  to  para¬ 
meter  variations  which  simulated  modal  error,  was  made.  Sys¬ 
tems  with  both  stable  and  unstable  controllers  were  considered. 


Investigation 

This  chapter  consists  of  discussions  of  the  results 
obtained  in  investigating  the  areas  presented  in  the  outline. 

Singular  Perturbation 
Optimal  Control 

This  section  of  the  investigation  applies  singular  per¬ 
turbation  techniques  to  reduce  control  spillover.  Initially, 
a  single  co-located  sensor/actuator  pair  were  used  to  control 
the  first  three  modes  and  suppress  the  fourth.  The  weighting 
matrix,  Qs,  was  increased  from  10°  to  1014,  and  system  response 
was  examined.  As  Q  approached  1014,  the  real  part  of  the 
eigenvalue  associated  with  the  suppressed  mode,  which  was  in- 
itally  positive,  approached  zero,  indicating  the  method  was 
reducing  control  spillover.  However,  the  real  part  of  the 
eigenvalues  associated  with  the  three  controlled  modes,  which 
were  all  initially  negative,  also  approached  zero,  indicating 
a  general  loss  of  control  in  the  system.  Since  the  method  did 
not  appear  suited  to  single  sensor/actuator  operation,  two 
sensors  and  actuators  were  then  used.  In  this  case  control 
could  be  maintained,  while  control  spillover  was  reduced, 
however,  as  can  be  seen  in  Figure  3,  little  suppression  of  the 
fourth  mode  is  realized  until  Qs  gets  very  large  (-  1011)  as 
compared  to  R  and  F  (=  1) .  After  this  point,  the  curve  flattens. 
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and  an  infinite  weighting  would  theoretically  be  required  to 
achieve  complete  suppression.  The  suppressed  mode,  in  this 
case,  had  nn  eigenvalue  of  1.501,  ±jl20.37  using  a  Q  of  1, 
and  an  eigenvalue  of  0.0021,  ±jl20.91  using  a  Qg  of  1014,  which 
demonstrates  the  reduction  in  control  spillover  achieved  using 
this  method.  It  can  be  concluded,  therefore,  that  while  this 
method  is  effective  in  reducing  control  spillover,  it  will 
never  actually  eliminate  it.  The  control  spillover  will  be  re¬ 
duced,  as  Qg  increases,  only  to  the  point  where  the  numbers  in¬ 
volved  become  so  large  that  computer  roundoff  error  becomes 
significant.  The  method  also  appears  to  require  a  reduction 
in  control  to  reduce  control  spillover,  explaining  why  control 
was  lost  in  single  actuator  operation.  These  results  can  be 
substantiated  in  light  of  Coradetti's  Limit  Method,  where  com¬ 
puter  roundoff  error  for  very  large  Qs  eventually  caused  prob¬ 
lems  in  matrix  inversion,  necessary  in  solving  the  matrix 
Ricatti  equation. 

The  Transformation  Method 

In  this  section,  a  transformation  matrix  is  generated  to 

i 

produce  a  gain  matrix  for  which  control  spillover  is  eliminated. 
Two  co-located  sensor/actuator  pairs  were  used  to  control  the 
first  two  modes  and  suppress  modes  three  and  four.  In  apply¬ 
ing  the  Transformation  Method,  Bg  should  be  made  not  full  rank. 
In  the  case  above,  the  Bg  matrix  has  two  non-zero  singular 
values  usually  associated  with  it,  determined  from  singular 
value  decomposition.  Since  the  number  of  non-zero  singular 

I 

■  values  of  a  matrix  equals  its  rank,  if  one  of  these  can  be  made 
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approximately  zero,  through  appropriate  actuator  placement, 
the  rank  of  Bg  would  decrease  from  two  to  one.  The  Singu¬ 
larity  Computer  Program  (Appendix  B)  was  used  to  produce  Fig¬ 
ure  4,  from  which  actuator  locations  could  be  determined  to 
decrease  the  rank  of  B  .  In  this  figure,  actuator  1  position 
is  held  constant  and  the  minimum  singular  value  of  Bg  is 
plotted  as  actuator  2  position  is  changed.  Actuator  combina¬ 
tions  are  chosen  such  that  the  minimum  singular  value  approaches 
zero,  excluding  co-located  actuators,  which  produce  trivial  sol¬ 
utions.  Two  cases  are  examined  in  this  figure  yielding  useful 
actuator  combinations  of  :  0.6,  0.285;  0.6,  0.89;  0.2,  0.585; 

0.2,  0.855.  The  method  appeared  to  work  extremely  well  using 
these  actuator  combinations,  essentially  yielding  zero  eigen¬ 
values  for  the  two  suppressed  modes  and  uncoupling  system 
eigenvalues  to  six  decimal  places.  Table  I  demonstrates  both 
of  these  results.  The  largest  control  spillover  terms  were 
on  the  order  of  10"4.  These  terms  would  decrease  as  actuator 
placement  approached  more  precisely  the  point  producing  a  zero 
singular  value.  A  comparison  between  a  system  utilizing  the 
Transformation  Method,  and  the  same  system  not  using  it,  is 
presented  in  Table  II,  demonstrating  the  Transformation  Method's 
effectiveness  in  suppression  and  uncoupling  of  system  eigen¬ 
values. 

Demand  on  the  system  was  then  increased  to  control  the 
first  three  modes  and  suppress  modes  four,  five,  and  six.  As 
pointed  out  in  the  presentation  of  the  Transformation  Method, 
actuator  placement,  when  there  are  more  modes  to  be  suppressed 
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Singular  Value  of  B  ,  2  Actuators,  2  Suppressed  Modes 


Table  I 


Modal  Suppression  and  Eigenvalue  Uncoupling 


Using 

the  Transformation 

Method 

Eigenvalues 

of  A  +  B  G 

-1.341328 

0. 

(C-l) 

— c  — c— 

-15.953688 

0. 

(C-l) 

-5.979912 

± j  20 . 493 

(C-2 ) 

Eigenvalues 

of  Ac  -  K  Cc 

-0.923323 

± j  3 . 591 

(0-1) 

-13.493944 

± jl5 . 910 

(0-2) 

Eigenvalues  of 

the  System 

-2 . 883098xl0~7 

± j 120 . 912 

(S-4) 

+3 . 058266x10“® 

±161.701 

( S— 3 ) 

-1.341328 

0. 

(C-l) 

-15.953691 

0. 

(C-l) 

-5.979912 

±j20 .493 

(C-2) 

-0.923323 

± j  3 . 591 

(0-1) 

-13.493942 

± jl5 . 910 

(0-2) 

Two  Controlled,  Two  Suppressed  Modes 

Sensor  Locations 

0.4 

0.2 

Actuator  Locations 

0.6 

0.285 

Singular  Values  of  B 

B  2'12 

0.000011 

(S-x)  -  (Suppressed-Mode)  Eigenvalues 
(C-x)  -  (Control led-Mode)  Eigenvalues 
(0-x)  -  (Observer-Mode)  Eigenvalues 


Uncontrolled  Structural  Damping  Ratio  =  0.0 
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Effectiveness  of  the  Transformation  Method 
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than  actuators,  could  not  be  chosen  such  that  one  of  the  non¬ 
zero  singular  values  of  Bg  would  go  to  zero.  However,  actuator 
placement  could  be  chosen  to  minimize  a  non-zero  singular 
value.  The  non-linear  behavior  of  this  singular  value  is  dem¬ 
onstrated  in  Figure  5  which  represents  three  of  the  actuator 
placement  combinations  used.  The  figure  was  generated  setting 
one  actuator  at  a  specific  location  and  varying  the  location  of 
the  second  actuator.  The  curves  represent  the  minimum  singular 
value  found  using  all  possible  actuator  locations  with  a  step 
size  of  0.05  from  0.0  to  1.0.  The  curves  also  demonstrate  that 
the  only  actuator  combinations  that  can  produce  a  zero  singu¬ 
lar  value,  indicating  a  reduction  in  the  rank  of  Bs,  is  where 
the  actuators  are  co-located,  representing  a  trivial  solution. 
The  minimum  value  of  these  curves  will  provide  the  best,  in 
the  least  squares  sense,  actuator  combination  to  be  used  with 
the  method.  Results  obtained  in  this  case  are  presented  in 
Tables  III  and  IV.  Table  III  represents  the  results  obtained 
using  the  actuator  combination  (0.95,  0.79),  which  produced 
the  smallest  singular  value  (0.16),  and  Table  IV,  an  actuator 
combination  of  (0.65,  0.90),  which  produced  a  singular  value 
of  0.265.  The  third  actuator  combination  (0.90,  0.66)  pro¬ 
duced  results  very  similar  to  those  in  Table  IV  and  was  not 

presented.  The  spectral  radius  of  B  G  is  the  square  root  of 

— s— 

the  sum  of  the  squares  of  the  elements  in  BgG.  Both  systems 
shown  did  very  well  in  controlling  and  suppressing  the  modes 
required.  Also,  uncoupling  of  system  eigenvalues  was  very 
good.  It  can  be  seen  when  comparing  these  two  tables  that 
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Singular  Value  of  B  ,  2  Actuators,  3  Suppressed  Modes 


Table  III 


System  Response  Using  the  Transformation  Method 
With  a  Minimum  Singular  Value  for  Bg  of  0.16 


Eigenvalues 
of  he  +  ®cG 


Eigenvalues 

of  A  -  K  C 
— c  —  — c 


Eigenvalues 
of  the  System 


-1.2101 

± j  3 . 334 

(C-l) 

-0.5948 

± j  22 . 022 

(C-2) 

-0.3755 

± j  61 . 701 

(C-3) 

-6.8487 

0. 

(0-1) 

-2.3939 

0. 

(0-1) 

-16.4721 

± jl5 . 624 

(0-2) 

-31.5344 

± j  39 . 903 

(0-3) 

-1.4926 

±j298. 525 

(S-6) 

-0.9993 

± jl99 . 852 

(S-5) 

-0.6039 

± j 120 .909 

(S-4 ) 

-1.2101 

± j  3 . 334 

(C-l) 

-0.5948 

± j  22 . 022 

(C-2) 

-0.3755 

± j  61 . 201 

(C-3) 

-6.8494 

0. 

(0-1) 

-2.3937 

0. 

(0-1) 

-16.4707 

± jl5 . 625 

(0-2) 

-31.5413 

± j  39.905 

(0-3) 

Three  Controlled,  Three  Suppressed  Modes 

Sensors/Actuators  Co-located  at  0.95  0.79 

Uncontrolled  Structural  Damping  Ratio  =  0.05 

(C-x)  -  (Controlled-Mode)  Eigenvalues 
(0-x)  -  (Observer-Mode)  Eigenvalues 
(S-x)  -  (Suppressed-Mode)  Eigenvalues 

Control  Spillover,  BgG,  Largest  Term  =  0.149 
Spectral  Radius  of  BgG  =  0.280 
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Table  IV 


System  Response 

Using  the  Transformation  Method 

With  a  Minimum 

Singular  Value 

for  B  of  0.265 
— s 

Eigenvalues 

of  A_  +  B  G 
— c  — c— 

-1.056 

-0.159 

±  j  3 . 377 
±  j  22 . 033 

(C-l) 

(C-2) 

-0.330 

± j  61 . 699 

(C-3) 

Eigenvalues 

of  he  -  £  9c 

-3.173 

±jl. 365 

(0-1) 

-19.417 

±j  21 . 484 

(0-2) 

-25.596 

± j  38.826 

(0-3) 

Eigenvalues 
of  the  System 

-1.492 

± j298 . 525 

(S-6) 

-0.999 

± j 199 . 852 

(S-5) 

-0.602 

± j 120 . 910 

(S-4) 

-1.056 

±33.377 

(C-l) 

-0.159 

± j  22 . 033 

(C-2) 

-0.330 

± j  61 . 694 

(C-3) 

-3.173 

ljl.366 

(0-1) 

-19.416 

±j21.483 

(0-2) 

-25.597 

± j  38.828 

(0-3) 

Three  Controlled,  Three  Suppressed  Modes 

Sensors/Actuators  Co-located  at  0.65  0.90 

Uncontrolled  Structural  Damping  Ratio  =  0.05 

(C-x)  -  (Controlled-Mode)  Eigenvalues 
(O-x)  -  (Observer-Mode)  Eigenvalues 
(S-x)  -  (Suppressed-Mode)  Eigenvalues 

Control  Spillover,  B  G,  Largest  Term  =  0.223 
Spectral  Radius  of  B^G  =  0.314 
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control  spillover  is  less  and  the  system  responds  better  as 
the  singular  value  decreased.  These  effects  were  verified 
through  further  testing. 

Observation  Spillover  Reduction 

This  section  uses  the  Transformation  Method  to  reduce 
observation  spillover.  Since  the  Transformation  Method 
apparently  worked  so  well  in  minimizing  control  spillover, 
the  method  was  applied  to  reduce  observation  spillover  with 
the  thought  that  system  response  may  be  further  improved.  In 
the  cases  tested,  the  method  was  applied  to  reduce  both  ob¬ 
servation  and  control  spillover  simultaneously.  Again,  two 
co-located  sensor/actuator  pairs  were  used  to  control  three 
modes  and  suppress  three.  The  procedure  was  identical  to 
that  used  in  reducing  control  spillover,  and  was  simplified 
due  to  the  fact  that  the  sensors  and  actuators  were  co-located. 
The  results  obtained  proved  interesting  in  that,  although 
observation  spillover  was  decreased,  the  method  actually 
appeared  to  decrease  the  overall  stability  of  the  system.  In 
the  cases  tested,  all  system  eigenvalues  remained  essentially 
the  same,  except  for  three.  These  three  were  observer  eigen¬ 
values  associated  with  the  ^  -  K  Cj  matrix.  A  representative 
example  illustrating  this  effect  is  shown  in  Table  V,  where 
control  spillover  is  reduced  in  both  cases,  but  observation 
spillover  in  only  one.  The  eigenvalues  of  the  three  suppress¬ 
ed  modes,  which  are  directly  influenced  by  the  combined  effect 
of  observation  and  control  spillover,  are  not  significantly 
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'changed,  however,  the  eigenvalues  of  the  Ap  -  K  Cc  matrix,  an 
indication  of  the  rate  at  which  the  observer  error  terms  decay, 
have  generally  shifted  to  the  right.  Evidently,  the  Transfor¬ 
mation  Method,  applied  in  reducing  observation  spillover,  re¬ 
duces  observability,  much  as  it  reduces  control,  in  achieving 
a  reduction  in  spillover.  This  reduction  in  observability 
has  the  effect  of  producing  a  gain  matrix,  K,  that  moves  the 
eigenvalues  of  Ac  -  K  Cc  closer  to  the  unstable  region.  From 
these  results  it  appears  that  applying  this  method  to  reduce 
the  observation  spillover  in  a  system  would  degrade  system 
response. 

Sensitivity  of  System  to 
Parameter  Variation 

This  final  portion  of  the  investigation  introduces  mode 
shape  errors  to  determine  system  sensitivity.  Systems  with 
a  stable  versus  an  unstable  controller  are  examined.  The 
system  considered  used  two  sensors  and  two  actuators  to  con¬ 
trol  the  first  two  modes  and  suppress  modes  three  and  four. 

The  Transformation  Method  was  applied  and  errors  in  modal 
amplitude,  evaluated  at  the  actuator  locations,  were  simulated 
up  to  ±20%  in  5%  increments.  Although  a  great  deal  of  data 
was  generated  for  systems  with  either  a  stable  or  an  unstable 
controller,  no  definite  determination  could  be  made  that  an 
unstable  controller  produced  a  system  that  was  more  suscept- 
able  to  modal  errors.  It  was  noted,  however,  that  as  the 
simulated  error  increased,  the  uncoupling  of  system  eigen- 
-  values  degraded  significantly  more  in  a  system  with  an  unstable 


controller,  than  in  one  with  a  stable  one.  Table  VI  repre¬ 
sents  the  results  typically  obtained  in  this  analysis.  The 
effect  of  a  20%  simulated  error  for  a  system  with  a  stable 
controller  is  compared  to  a  system  using  an  unstable  controller. 
The  same  analysis  was  conducted  on  a  system  controlling  the 
first  three  modes  and  suppressing  modes  four,  five,  and  six, 
producing  similar  results.  Although  the  eigenvalues  of  the 
system  did  not  vary  as  significantly  as  the  uncoupling  effect 
in  this  model,  it  is  possible  that  this  uncoupling  sensitivity 
could  adversely  affect  a  system  representing  a  more  complex 


model . 


Uncoupling  Differences  for  a  Stable  and  Unstable  Controller 
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Conclusions 


The  two  main  conclusions  to  be  drawn  from  this  report 
are  that  the  destabilizing  effect,  caused  through  interaction 
of  observation  and  control  spillover,  on  suppressed  mode  re¬ 
sponse  can  effectively  be  reduced,  and  that  the  system  examin¬ 
ed  was  quite  insensitive  to  modal  error. 

Singular  Perturbation  Optimal  Control  did  not  yield  re¬ 
sults  as  encouraging  as  those  found  using  the  Transformation 
Method,  however,  the  method  was  effective  in  reducing  control 
spillover.  Problems  arise  in  matrix  inversion,  due  to  com¬ 
puter  roundoff  error,  when  weightings  used  become  very  large, 
and  an  infinite  weighting  would  theoretically  be  required  to 
completely  eliminate  control  spillover  effects.  Only  in  this 
limiting  case  does  the  Singular  Perturbation  Method  yield 
results  comparable  to  those  obtained  using  the  Transformation 
Method  with  appropriate  placement  of  actuators. 

The  Transformation  Method  was  found  to  be  very  effective 
in  eliminating  control  spillover  and  uncoupling  of  system 
eigenvalues,  when  the  number  of  actuators  in  the  system  was 
equal  to/or  greater  than  the  number  of  suppressed  modes. 

Where  the  number  of  modes  to  be  suppr-  ssed  exceeded  the  number 
of  actuators  in  the  system,  the  method  provided  the  best,  in 
a  least  squares  sense,  actuator  locations  to  minimize  control 
spillover. 
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Reduction  of  observation  spillover  using  the  Transfor¬ 
mation  Method  was  found  to  degrade  system  response.  The  re¬ 
duction  of  observability,  an  adverse  effect  of  decreasing  the 
observation  spillover,  apparantly  was  sufficient  to  alter 
the  observation  gain  matrix,  such  that  the  eigenvalues  of 
Ac  -  K  Cc ,  which  indicate  the  rate  of  decay  in  observation 
error,  shifted  to  the  right,  producing  a  destabilizing  effect 
on  system  response. 

The  modal  error  analysis  conducted  did  not  substantiate 
a  presupposition  that  a  system  with  an  unstable  controller 
would  be  more  sensitive  to  parameter  variation  than  one  with 
a  stable  one.  In  both  cases,  it  was  found  that  the  system 
was  quite  insensitive  to  modal  error,  with  the  eigenvalues 
of  the  system  never  differing  greater  than  the  errors  in¬ 
duced.  Data  generated,  however,  did  show  that  the  uncoupling 
of  system  eigenvalues  was  significantly  more  sensitive  to 
parameter  variation  in  a  system  with  an  unstable  controller. 
Although  this  uncoupling  sensitivity  did  not  affect  system 
response  in  the  model  tested,  it  remains  to  be  determined  if 
it  would  cause  problems  in  a  more  complex  model. 
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Recommendations 

In  a  large  complex  structure  the  number  of  modes  to  be 
controlled  and  suppressed  could  become  quite  large.  As  the 
number  of  modes  considered  in  a  model  increases,  the  computer 
capabilities  as  well  as  sensor  and  actuator  requirements  to 
control  the  system  would  have  to  increase.  The  first  concern 
therefore,  in  modelling  a  structure,  is  which  modes  should 
be  considered  critical  and  included  in  the  system.  A  method 
of  determining  these  critical  modes,  in  a  complex  structure, 
would  have  to  be  the  first  step  in  implementing  the  control 
methods.  Further  study,  of  the  methods  presented  in  this 
thesis,  on  a  complex  model  is  needed  to  determine  the  viabili¬ 
ty  of  this  approach  to  controlling  large  space  structures. 
Unstable  controllers  in  a  complex  model  should  also  be  ex¬ 
amined  for  sensitivity  to  modelling  inaccuracies,  as  the 
analysis  presented  could  not  dismiss  this  possibility. 
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Appendix  A 

Main  Computer  Program  Listing 
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PROGF AH  lUOBS2(  INFUT=/Fu, OUTPUT*/ 132, TAPE5 , TAPES, TAPE?) 

DIMENSION  3l  (  2  u  )  i  G  (  2  >  2'  )  ,  C  (  2 »  2  u )  »B(j|  2)  ,  W  >i  ft  T  (20)  »Z(6»s)  »W<n3) 
•/.,AFk1ME(23,2u)  ,T(*/J,*0)  ,W(5G)  ,AP(f0,:50),wr  ((•,&>  ,X  (6,6)  ,S  (6,6)  , 
*/.A<(2i  ,2)  ,Za(2),  2C )  ,Ze(23 ,2i>  ,  AD  (2C,2L  >  ,2011(20,20  ,C012(20,2D>  # 
*/.CK(2,  1C)  ,BR(b,lI>  ,AR(1<-  ,13)  ,H12(>C,10)  ,H  21  (10,20  ,ZZ(23,20)  , 

*4  AO  AMP  (20,2))  ,  ADA-PR(  20, 20),  A  (6,  3)  ,  3RSIME,  2)  ,BIN(6,b)  ,  KT(2,b)  , 
‘/.SIMA  )  ,T1  (2,1)  ,R(  2,2),rIN(2,2)  ,Gl  (1,6)  ,*09IN<2,2)  ,  ROB  (2, 2)  ,31  (2,6) 
*/.,Cl  «  ,2)  ,CTR(5,  2)  ,CTRC(6,6)  , F03(>,9)  ,  XOB  (  5 , 6  >  ,  ZOB  ( b ,  6)  ,RC(2,6)  , 

XRT (1, 1),T1T(1,2),71TR(1,2),B1T(6, 1), *T BT ( L ,6) ,RT INV( 1, 1) ,31TT (1,5) 
X,Z5(fc,6),8X(6,2),BRX(6,2),CX(2,6»  ,  CSX  (2,0  ,  BRISK  6 , 2)  ,  T2  ( 2, 1) 
*/.,CIN(fc,6)  ,CkSIN  (:  ,2)  ,STN2(4)  ,  T 2T3  ( 1 , 6)  ,  KTi  (1 ,6)  ,  CRT  (6 , 2)  ,T2T(1,2) 
CONMCN/SHAPE/u,  M 
COhMG N/FUNCT/4P IIHE,  NF 
COHMUN/HAINl/.(DIM,fniHl,COHl 
COMMC  -N/INOU/KIN ,KCJT ^UNCH 
COMNON/MA IN2/C0H? 

EXTEF  NAL  FN 
kEAL  L,M 
REAL  KT 
COHPLEX  W,Q 

DATA  IB/c/, IS/u/,  IX/6/,lG/2/,IZ/5/,I&DP/2t /, xO/S 0/ ,IZZ/20/, 

•/.IZA/2  0/,  I9R/5/,  H  21/  1C/ ,  i A  P/5  C/ , T  AK/ 2  it /,  1 2 /2/,  ICR/2/,  I H12/20/ , 
*/.IZE/2  0/,IAR/13/  ,IARIKV/1C/,1AIBR/10/ 

NC  =  i 
NR  =  3 
N=  3 

NC2=2*NC 
NR2=NR"»2 
NT  =NR  *N 
N2=NCIH=2*N 
NF  =  2*  N2+NR2 
NF  0=2*N2 
N01N1=N0IH+1 
ZETA=C.0u5 
IN1T*  0 

1  REAO*  ,  L,  M, PAI  Ri,  P  AIR?,  «=AIR3,  PAIR*,  OKI)  ,I=1,N')  ,EI 
IF  (EG  F  (5L INPUT)  •  N  E,  9  • )  STOP  **GRAI  EFJL  L  Y** 

READ*,  (WT(I,I)  ,I=l,N2) 

DO  b  1*1, N2 
00  5  J*l, N2 
FO 6 (i  ,J)=0. 

6  CONTI  NUE 

READ*  ,  (FOB(I,I> ,1=1, N2) 

R(  1,1 )  *R(  2,2)  =ROR  (1,  1)  =F;09(  2, 2)  =1  • 

Ml,2)=R(2,l)=R0R(l,  2)=R09(2,1)=). 

00  3  1*1.  NT 

3  WNAT  ( I  )  =  (  (BL  (I )  /L )  2)  *■'  SORT  (  El/Ml 

00  H  1*1, N2 
00  4  J=1.N2 
A(1,J)=0. 

AO  (I,  J)*u  • 

AO AMP (I, J)=C. 

AR  (I.  J)  =  0. 

AOaHPRU,  J)*0. 

CONTINUE 

A  (4 , 1 )  =A0  (4,0*  A3  AMP  (4 , 1)  =-WNAT  (1)  **2 
A(t?,2)=A9(5,2)=ADAHP  (5,2)  =-WNAT  (2)  **2 
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ACt,2)=Ac  (•d,3)  =  A:AHP  (6,3)=-WNAT  (3) **2 
A  ( 1  »*•  )=AD(1,‘4)=A’)AMP  (i,4)=l. 

(  A(2,5)=AD<2,‘ >  =  AnAM^>  <2,5>=1. 

A<3,c)=A0(3,'.)=A:>A4»  <3,6)=1. 

A  ( •*,<•  )  =AD  (4,  4 )  =  ADAMP  (4,4)=-2.*7Er£*HsAT  (1) 
A(i,f  )=AD(5»F)=ADA.'10  (5 , 5  )  =-2  .  *  Z  ET  A*  UN  AT  (  2  ) 
A(t,C  )=AO(i,c)  =  AHA'IP  (6,6)=-2.*2Er^WNAT  (3) 
Ar<  ( 4 1  i)=A0AMPR(4,  i)  =  -WNATK)  **2 
AR<5,  2>=ADAMPR(5,  2)  =  -WNAT  (5 )  **  2 
AR<6,  3)=ADAMPR(S,  3)  = -WMAT  (o)  **  2 
AR  ( 1,  h  )  =A0AMPR  (  1,  *t )  =  1. 

AR  (2,  5) = AOAMPR (  2,5)  =  1. 

Ak(3,6)=ADAMP^(3,£>  =  1. 

AR  (4,  h)  =  AOAMPR(4,  ‘  >  =  -2.*ZETA*WNU  (4) 
AR<t»,5)=A0ArtPR(3,’  )  = -2. * ZET A* WNAT  (3) 
AR(o,6)=ADAMPR(S,t>)  =  “2**Z£TA*WNAr  (5) 

PRINT  liij,  FAIRl»PAIR2,PAlR3,PAIR+jL,M,EI 
PRINT  ill  5,  <3L<l),I=i,NT),<WNAT(I),I=l,NT) 
00  V:  1=1,  N2 
B(i,l )=3(I,2) =  u  . 

13  CONTINUE 

DO  12  1=1, N2 
00  12  J=i,N2 
IF (I . Nl. J)  WT (I, J)  =  3  , 

12  CONTINUE 
J=N+1 

DO  15  I=J,N2 

.  CALL  MODESCBCI, 1) ,BL (I-N) ,PAIRi> 

(  CALL  MC0ES(9(I,2) ,9L (I-N) ,PAIR2) 

15  CONTINUE 

DO  35  1= J  ,N2 
C(1,I)=C<2,I)=0  . 

35  CONTINUE 

00  3 #  1=1, NR2 

CR  (I,  I)  =BR(1 , 1)  =CR(2,I)=BR(I,2)=3. 

3?  CONTINUE 

DO  4  u  1=1, NR 

CALL  MOO£S<C(1,I> ,Bl (I) ,PA1R3) 

CALL  M0DES(C(2,I)  ,3L  (I) ,PAi*4> 

4  .■  Continue 

00  4f  1=1, NR 

CALL  MOOES (BR (i  ,  1)  ,BL  CI+N)  ,PAI?1> 

CALL  MOOES (8R Cl  +  N,2)  ,BL<I*N>  ,PAI*2) 

CALL  MC0£S(CR(1,I)  ,3L(T+N)  ,PAlR-*> 

CALL  MC0ESCCRC2 ,1) ,3 ICI+N)  ,PAIR4> 

45  CONTINUE 

00  53  1=1 ,NR2 
00  53  J=i,NC2 
CRSIN (I, J )=CR( J,  I ) 

BRSIN  (I , J ) =BR (X , J ) 

53  CONTINUE 

00  &«•  1=1,  NR2 
00  54  J=1,NR2 
(  CINCI  ,  J)  =  0* 

BIMCI  ,  J)  =  0. 

54  CONTINUE 

00  50  1=1, NR2 
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CINC  ,1>  =  1. 

BIN(1  ,l)  =  i. 

66  CONTINUE 
PRINT  16-*  5 
PRINT ♦  ,  "  F  " 

00  tl  1=1, N2 

PRINT  *  ,  (WT  (I*  i)  ,J=1,N2) 

61  CONTINUE 
WRITE*,  "  " 

PRINT*,  "  P.  M 
00  fc2  1=1,  NC  2 

62  Pi\lNT*,(R(I»J)  ,J=i,NC2) 

WRITE*,  M  " 

PRINT*,  ’*  FOB“ 

00  63  1  =  1, N2 

63  PRINT*,  (  FOB  (I ,  J)  ,J=1,N2) 

WRITE*,  "  ** 

PRINT*,  "ROB" 

DO  6«-  1=1,  NC2 

61*  PRINT*, (R0B(I,J),J=1,NC2) 

PRINT  12 JG 
DO  *5  1  =  1  ,N2 

PRINT  1 3c-  0 ,  (A DA MP  (I ,  J),J  =  1,N2) 
i-9  CONTINUE 
PRINT  1325 
DO  51  1=1, NR2 

PRINT  135  J,  (AOAHPRd  ,J)  ,  J=1,NR2) 

'  CONTINUE 
PRINT  141/0 
DO  5E  1=1, N2 

PRINT  1330,C3(I,J)  ,  J=1,NC2) 

PRINT  l<+y  0 
DO  59  1=1, NR2 

59  PRINT  135G,  (BR(I,  J>  ,  0=1, NC2) 

PRINT  15  JO 

PRINT  13wO, (C(1,I),I=1,N2) 

Pi\INT  13jO,(C(2,I),I  =1.,  N2) 

PRINT  155  0 

PRINT  13^0, <CR (1,1) ,1=1, NR2) 

PRINT  13v  3,  <CR(2,I)  ,1=1,  NR2) 

PRINT  15 1*5 

CALL  LSV0FC9R3I N, NR2 ,NP2 ,NC2 , BIN, NR2,  NR2*  SIN ,WK, IER) 
DO  6P  1=1  ,NC2 
WRITE*,  SIN(I) 

68  CONTINUE 

00  121  1=1,  NC2 
Tl (I, 1)=BRSIN(I ,2) 

121  TIT (1 ,I)=Tl(I,l ) 

PRINT  156  0 
00  67  1=1 *NC2 

67  WRITE*, TlCI.l) 

00  97  I*1»NR2 
00  97  J*1 »NC2 
CT  (I,  J)=C  (J, I) 

97  CONTINUE 

CALL  LINV2F<RJ8, NC2, NC 2, ROBIN,! ,<K,IER) 

CALL  VNULFF(CT,  ROBIN ,N2 , NC2, NC2, i 2, S3 2 ,CTR,N2,IER> 
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CALL  VMULFF(CTR,::,N2|NC2,N2,N2,N;2,Cf  RC,N2,IER) 

CALL  MfcIC  (N2,AtC7r.C,  FO^jXOB,  ZOB,  i  R,  I  V  IT ) 

CALL  VMULFF <R33IM,C,  NC2  ,  NC2,U2,  NC2,  TC  2,  RC,  NC2,I  ER> 
CALL  VMULFF(RC,XGB,NC2,N2,N2,NC2,  N2, <T,NC2 ,IER) 
WRITE*#  "  M 
PRINT  23*  0 
DO  9f  I=1#NC2 

9S  PRINT  13-0#  < < T  (I»J)  #J  =  1#N2) 

00  122  1=1,  NR2 
00  12  2  J=1,NR2 
S ( 1 j J )=u . 

122  CONTINUE 

CALL  VMULFF  (  B#T  1,  NR2  , NF  2, NC,  NR2, VC2,  5 IT,  NR 2,  IEf< ) 

DO  111  1=1, NR2 
BlTT( 1,1) =91 T  <1,1 ) 

161  CONTINUE 

CALL  VHULFF(T1T  ,NC  ,NC2  , NC2 ,  NC , 'J C2,  T  1TR , NC,  1E«.) 

CALL  VMULFF(T1TR,T1,NC,NC2,NC ,NC, KC2,RT »NC»IER) 

CALL  L1NV  2F (RT ,  NC,NC  ,RT I NV, 0 , WK , I ER> 

CALL  VMULFF  (BIT  ,  RTIN  V,  Nk2, NC ,  NC  ,  VF.2,  SC ,  9KI N»  NR2,  IER) 
CALC  VHUlFFTBRIN, BIT T, NR2 ,NC , Nk2, NR2 , NC , $ , NK2, I ER) 
INIT=  0 

CALL  MKIC (NR2,A,S,WT ,X,Z,MR, INITJ 

IF  (MR  •  EQ.  ”1)  STOP  "MRIC  DID  NOT  CONVERGE" 

CA LL  VMUL  FF ( RTI NV , 31 TT , NC  ,NC , NR2,  NC , VC , RT  PT , NC, I ER) 
CALL  VMULFF  (RTBT,X,NC,NR2,NR2,NC, NR2, Gl,NC,IER) 

IF  TIER. EQ. 129)  STOP  "BT*X  BAD" 

CALL  VHULFF(T1,G1,NC2,  NC,NR2,NC2,  NC ,  3  ,  NC  2 ,  IER) 

PRINT  16C 0 

PRINT  13Ju,(G(l,I) ,I=1,N2) 

PRINT  130  3, (G(2,I) , I =1, N  2) 

CALL  VMULFF (P,G,N2,NC2jN2,IB,IG,Z#IZ#IER) 

IF  (IER. EO.  129)  STOP  "B”  G  BAD** 

DO  it  2  1=1,  N2 
DO  1C  2  J=1,NC2 
IE  2  AK (I, J)  =  KT ( J, I) 

CALL  VMULFF <AK,C,N2, NC2 , N2, IAK,IC , Z3 , IZ9 , I ER) 

IJOB=  C 

CALL  E1GkF(A0AH»,N2,  IAPP ,  I  JO  B  ,  W,  3 , 10,  WK,  IER) 
IFCIER.NE.O)  PRINT* , "TER  =  ",IER 
PRINT  90b 
00  57  1=1, N2 
PRINT*, W(I) 

57  CONTINUE 

DO  6fc  1=1, N2 
DO  66  J=i  »N2 
6b  CONTINUE 
PRINT  175  0 
00  <JC  1=1, N2 

80  PRINT  130  0,  (Z5(  I,  J)  ,  J=1,N2) 

I JOB= 0 

CALL  E1GrF(Z5,N2, N2,IJCB,W,Q,IQ,(X,IER) 

PRINT  1765 
00  81  1=1, N2 
PRINT*,  W(I) 

81  CONTINUE 
PRINT  1760 
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FACTC Rl=  3 ,5 
F4LTC'R2=1.2 
WRITE*,  »  " 

WRITE*  ,  “  ** 

WRITE*,  M  " 

WRITE*,  ”  “ 

WRITE*,  M  " 

CO  120  I J=i,9 

PRINT  * ,  ••FACT  OR  1  =  ”,  FACT OR1,  ”  FACTOR  2  = 

00  11 M  1=1, N2 
00  111,  J=1,NC2 

3RXC1  ,J)  =  3X(I,J)  =  CX<  J,I)=CRX(J,I>  =  d. 

11*  CONTINUE 

BX  (£, i)=FACT0R2*8 <5, 1) 

BX  <4,  1)  =F ACT0R1*9  (4,  1) 

BX(6,l)=FACT0Ri*r(6,  1) 
ex  (6,  2)=FaCTGR2*3(5,  2) 

BX (5, 2)=FACT0R1*3 (3, 2) 

BX  (*  ,  2)=FACT0*2*P  <4,  2) 

?RX<?  ,  1)  =  FACT 0R2*  3R(",1) 

BRX  (E  ,  2)  =  FACT  ORi  *  BR<9,2) 

BRX  (C  ,  1)  =FACT OR  1*  BR(6,i) 

BRX(s  ,2)  =FACTJR2-BR(6,?) 

PRX  (~  ,  1)  =  FACT  OKI*  HR(a,1) 

BRX  (»•  ,  2)  =  FACT  0R2r  3R  (  *,  2) 

CX  <1, 1>=FACT0R2*C Cl, 1) 

CX  (1,  2)=FACTGU*C  (1,2) 

CX (2, 1 ) =F  ACT0R1 * C (2, 1) 

CX (2, 2)=FACT0R2*C (2, 2) 

CX (1, 3)=F ACT0R2*C  (1,  3) 

CX (2,  3)=FACTGR1*C (2, 3) 

CRX  (1  ,l)=FACT0R2rCR(i,l) 

CRX<1 ,2)  =  FACT0R1‘CR(  1,2) 

CRX  (2 , 1)  =  FACT 0R1*  CR(  2,1) 

CRX(2,2)  =  FACTDR2'  CR(  2,2) 

CRX(1,3)=FACT3R2‘  CR(1,3) 

CRX  (2,3)=FACT0Rl'iCR(2,  3) 

CALL  VMULFF( PRX ,G ,N 2 ,NC 2 , N2, IBR,[G,H2l,IH21,lER) 
IF (1ER.GT  .129)  STOP  "Br*G  3A0" 

PRINT  I7i  0 
DO  111  1=1, N2 

PRINT  13w  0,  (H21  (I  ,  J)  ,  J=l, N2) 

111  CONTINUE 

CALC  VMULFF (BX,G,N2,NC2,N2,IB,IG,  7,IZ,IER) 

DO  5;  1=1, N2 
DO  SI  J=1  ,N2 
APRIhEd,  J*N2)  =  -7(I,  J) 

ZA(I,J)=A(I,J)-7(I,J) 

Z( 1 , J ) =A0  (I»J)-7(I,J) 

APRIP  E(I,J)=Z(I,J) 

ZZ (I, J)=Z  (I, J) 

APRIME(I  +  N2,J)  =  3, 

55  CONTINUE 
IJOB=  0 

CALL  E1GRF(ZZ»N  2, 177  ,1  JOB,W,  Q,IQ,  WK,  I  ER) 
IF(IER,NE.J)  PRINT*, ”IER  =  ”,IER 
PRINT  17 1 0 


,  FACT0R2 
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(00  6u  1=1, N2 

PRINT  131 0«CZCI,J),  J=1»N2) 
r  \  CONTINUE 
V  PRINT  180 C 

00  hi  1=1, N2 
PRINT  *,WC I) 

65  CONTINUE 

CALL  VMULFFUK,  CX , N2  ,NC 2 ,  N2, 1 AK ,  t  3 , Z3 ,179 ,  IEP.T 
00  1(  0  1  =  1, N2 
00  ICC  J*1,M2 
Z9CI,J)=A0(ItJ>-79(I,J) 

2A  (I,  J)  =  Z8(I,J) 

APKIkE  CI  +  N2,  J+N2)  =  Z9  (I,  J) 

1C',  CONTINUE 

uoa=o 

CALL  EIGRFCZA, N2* I7A ,1 JOB ,W,  Q, 19,  WK,  1 Ek) 
IF(Itk.Nt.O)  PRINT*-,  “If  R  =  “ ,  IER 
PRINT  2210 
00  U 5  1=1, N2 

PRINT  13L0»(ZB(I,J),J=1,N2) 

ICE  CONTINUE 

PRINT  23C  0 
00  110  1=1,  N2 
PRINT  * , W( I) 

11C  CONTINUE 

CALL  VPULFF ( A<,  CRX,  M  2,NC2,N2,  IAK,  ICR,  H12,1H12,IEK.) 
DO  112  1=1,  NR2 
00  112  J=1»N2 

f  APRINE(I+12,J)=H21CI,J) 

^  APKIME  (1*12,  J*b>  =  H21  (I ,  J) 

APRIht  (J,  1+12)  =  J. 

APRINE (J +6,1+12)= HI 2 CJ,I) 

112  CONTINUE 

00  113  1  =  13,  NF 
DO  113  J=13,NF 
APKINECI,  J)=ARCI-12,  J-12) 

113  CONTINUE 

00  115>  1=1, NF 
00  lib  J=1,NF 
AP (I, J)=APRIM£ (I, J) 

1U  CONTINUE 

PRINT  255  0 
00  95  1=1, NF 

99  PRINT  26j  0, CAPRINE (I ,J) , J=i, 6) 

WRITE* ,  "  M 
WRITE*,  M  •* 

00  126  1=1, NF 

126  PRINT  26JQ,  CAPRIHE(I,J)  ,  J=7,12> 

WRITE*,"  " 

WR*TE*,"  " 

00  127  1=1, NF 

127  PRINT  26CiO,(APRIME(I,J)*J=13,16) 

I JOB=  C 

/  CALL  EIGRF(AP,NF,IAF,IJOB|H,Q,nfWK,  IER) 

V  IFCieP.Nt.O)  PRINT*, "I^R  =  ",1ER 

WRITE*,  "  " 

WRITE*,  "  " 
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< 


9Cw 

111: 


12'.  0 
13.: 
1325 
135  d 
14  3  d 

14  yt 

15  ui? 
155», 
15  '.  v 
157  6 
1580 
15  SI 

15  0  0 

16  35 
lb^-3 

17  CG 
175  O’ 
17  $  s« 
1766 


1770 
18  CC 
1"C 

(  u 
22C0 

2300 


PRINT  24.0 
DO  1«.  5  1=1,  NF 
PRi.Nl  *  ,W(I) 

com:  nue 

FAGTl  K1=F  ACT0R1  ♦  .C  5 
FACUR2=FACTCR2-f;.f.5 
CONTINUE 
GO  TO  1 

FORMAT  (//10X, "THE  EIGENVALUES  OF  A "/» 

FOf  MAT(1H1/////  ,23X,  "CANTILEVER  3  EAH  WITH  TWO  SENSORS, 

7.2  ACTUATORS,  7  OCNTROLl  ED  MODES  AND  3  RESIDUAL  MODES", 

7//LIX  ,**Fj.RST  ACTUATOR  POSITION*", Fit . 8, 

7//60X , “SECOND  ACTUATOR  POSITION  =”,F14.8, 

7//5 L'X,"  FIRST  SENSOR  POSITION  *“,F1,.8, 

7//6G X,  “SECOND  SENSOR  POSITION  =*,P1V.3, 

7//5UX  ,“BEA.H  LENGTH=” ,F5 . 1, 

7/5GX,“8e.«M  MASS  DENS  ITY/UMIT  LENITH  =  “,^3.6, 

7/5  LX ,  “EE  A  M  STIFFNESS  (£1)  =  ”,F}.2) 

FORMAT  </i>  OX,  “RO  0TS  OF  "HE  f  HARACT  EP.IS  T  IC  EQUATION*”, 

7/>4X, “FIRST  MODE:  ",  E2r,,l*+,/55X, 'SECOND  MODE:  ", E20.lt, 

7/56X, “THIRD  MODE:  “,E2:.14, 

7./35X,  “FOURTH  MODE  :  ",E2i.l4, 

7/5  5X,  “FIFTH  MODE  :  “,£20.14, 

7/51X, “SIXTH  MODE  :", E2C.14, 

7/5 CX, “NATURAL  FREOUE NCI ES t“, 

7/55X, “FIRST  hJDE t  “, E2C .14, /55X, “SECOND  MODE:  “,E20.l4, 

7/5 5 X,  "THIRD  MODE:  ",E2:.14, 

7/55X, "FOURTH  MODE  :“,E20.14, 

7/ifcX, "FIFTH  MOO  E  :  “,E20.14, 

7/5EX,  “SIXTH  MODE:  “,E20.i4> 

FORMAT (//25X, *  THE  "A“  MATRIX*/) 

FORMAT  (/2X,6E16.;') 

FOF.MAT(//2dX,*THE  "A"  MATRIX  ( RE  SID  J  AL)  *  '  ) 

FOF.MAT  (/2X,6Elo  .7) 

FORMAT  <//2CX, 4  THE  “9“  MATRIX*/) 

FORMA  T(//15X,*T  HE  ”3"  MATRIX  FOR  RESIDUALS'  /) 

FORMAT  (//2JX,4THE  ”C“  MATRIX*/) 

FORMA  T(//21X,4T  HE  “C "  MATRIX  FOR  RESIDUALS*/) 

FOrvMAT  <//20X, 'SINGULARITIES  OF  THE  "5“  MATRIX  FOR  RESIDUALS*/) 
FORMA T(//2IX, *SINGUL  ARI1 I ES  OF  THE  C  MATRIX  FOk  RESIDUALS4/) 
FORMAT  <//2TX,*THE  “T  “  MATRIX  FOR  PR*/> 

FORMA T (//20X,4THE  "T“  MATRIX  FO?  CR*/) 

FORMAT  <25X,4THR  "G“  MATRIX4/) 

FORMAT  (//2LX,*THE  3R  TIMES  T  MATRIX4/) 

FORMA  T  (//10X,*THE  WEIGHTS  USED*/) 

FORMAT  <//57X,4T  HE  "A  *9G"  MATRIX*') 

FORMAT  (//23X, 'THE  A  ♦  EG  -  KC  MATRIX4/) 

FORMAT  <//lJX*4THE  EIGENVALUES  OF  A  +  BG  -  KC  MATRIX*/) 

FOrMA  T  (//1-JX,  "TO  ACCOUNT  FOR  MOOU  ERROR  THE  9  AND  C" 

7//1CX , "MaTRICIES  ARE  VARIED  BY  FiCTDR  1  FOR  ODD  ELEMENTS  AND  “ 
7//1CX , “3Y  FACTOR  2  F CR  EVEN  ELEMENTS"/) 

FORMA T(//10X,4BR  TIMES  G  MATRIX*') 

FORMA  T(//10X, "THE  EIGENVALUES  OF  A+3G"/)  . 

FORMA  T(10X, "THE  SMALLEST  REAL  PA?T  =  **,E23.13/) 
FQPMAT<//35X»*TH£  ”K”  MATRIX  TRAVS°35ED*/) 

FOrMA  T(//57X,*T HE  "A-KC"  MATRIX*') 

FORMAT  (//10X, "THE  EIGENVALUES  OF  A-C"/) 


FORMAT  (//57X,*T  HE 
FORMAT  <//23X, 'THE 
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24l-a  FOF.MA  T  (//1GX«“THE  eigenvalues  of  the  SYSTEM**/) 
25t-"  FORMAT  <//4uX,*‘  THE  SYSTEM  MATRIX*/) 

26  •'5  FORMAT  </2X,6E16. 7) 

END 
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SUBROUTINE  MODES<PHI,nLI,X> 

COhMCN/SHAPE/L*  M 
REAL  L,M 
BX  =  BL  I  /L'  X 

PHI  =  (  CSIrt(BLI>-SlNH(3LI>>M3IN<9<>-SINH(P<>IMCOS(3Ll>*COSHCBLl>>* 
y.(CLS(BX)-COSH(BX>  ))/  (SCRT  (M*  L)  *  SIN  (3LI)  *'  SXNH  (BLI )  > 

RETURN 

END 
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Appendix  B 

Singularity  Computer  Program  Listing 
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PROGRAM  LUOBS2  (INPUT  =  180,  OUTPUT  =/132,  TAPE  5,  TAPE  6,  TAPE  7) 
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16  CONTINUE 
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